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by 

Rodrick V. Chima 

NASA Glenn Research Center, Cleveland, OH 

and 

Timothy R. Conners and Thomas R. Way man 
Gulfstream Aerospace Corporation, Savannah, GA 

A computational analysis of a Gulfstream isentropic external compression supersonic inlet coupled to a 
Rolls-Royce fan was completed. The inlet was designed for a small, low sonic boom supersonic vehicle 
with a design cruise condition of M= 1.6 at 45,000 feet. The inlet design included an annular bypass duct 
that routed flow subsonically around an engine-mounted gearbox and diverted flow with high shock 
losses away from the fan tip. Two Reynolds-averaged Navier-Stokes codes were used for the analysis: an 
axisymmetric code called AVCS for the inlet and a 3-D code called SWIFT for the fan. The codes were 
coupled at a mixing plane boundary using a separate code for data exchange. The codes were used to 
determine the performance of the inlet / fan system at the design point and to predict the performance and 
operability of the system over the flight profile. At the design point the core inlet had a recovery of 96 
percent, and the fan operated near its peak efficiency and pressure ratio. A large hub radial distortion 
generated in the inlet was not eliminated by the fan and could pose a challenge for subsequent booster 
stages. The system operated stably at all points along the flight profile. Reduced stall margin was seen at 
low altitude and Mach number where flow separated on the interior lips of the cowl and bypass ducts. 
The coupled analysis gave consistent solutions at all points on the flight profile that would be difficult or 
impossible to predict by analysis of isolated components. 
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Objectives 

1. Develop CFD tools for coupled inlet / fan interaction. 

2. Apply the tools to a Gulfstream supersonic inlet coupled to the fan 
for a Rolls-Royce engine. 

3. Determine the performance of the system at the design point. 

4. Assess the operability of the system over the flight profile. 



Inlet / Fan for a Quiet Supersonic Aircraft 


R. Chima, NASA GRC, Sept. 2009 2 






AVCS code 



Axisymmetric Navier-Stokes duct analysis code by D. Tweedt and R. Chima 

• Node centered finite-difference formulation 
•AUSM + upwind differencing 

•Axisymmetric blockage term used to model gearbox blockage 

Explicit Runge-Kutta solver 

• Variable At 

• Implicit residual smoothing 
Turbulence model 

•Wilcox 2006 k-co turbulence model and stress limiter 
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SWIFT code 



3-D multiblock Navier-Stokes turbomachinery analysis code by R. Chima 

• Node centered finite-difference formulation 
•AUSM + upwind differencing 

Explicit Runge-Kutta solver 

• Variable At 

• Implicit residual smoothing 
Turbulence model 

•Wilcox 2006 k-cD turbulence model with stress limiter 

Coupled to AVCS with a steady mixing plane model 

•SWIFT blade-to-blade solution perturbed about AVCS solution using Giles’ 
characteristic BC 


Inlet / Fan for a Quiet Supersonic Aircraft 


R. Chima, NASA GRC, Sept. 2009 4 




Code Coupling - SYNCEX 





■ ■ ■ ■ ■ ■ ■ i 

halcyon: /home/fsrod/syncex 
File Edit View Terminal Ta_bs Help 

0 [7] syncdir > syncex & 

[1] 3429 

[8] syncdir > syncex: version 4.2 
syncex: accepting connections on port 58000 
[8] syncdir > 




$ 


halcyon: /home/fsrod/inlet 
File Edit View Terminal Tabs Help 


_ I 

pp ip 

Ur"lr I 

▼ 

-\ 


[5] inletdir > avcs < input > output & 


File 


▼ 

- ! - 

halcyon: /home/fsrod/fan 

Edit View Terminal Tabs Help 


i-tin x 


/ v 


[3] fandir > swift < input > output & 








/ 


SYNCEX - C language code written by Dr. Dan Tweedt of AP Solutions, Inc. 

• Runs in the background and handles data communication, storage, and 
synchronization between CFD codes 

• User routines read and write boundary condition data to SYNCEX 

• General interpolation routines provided 
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Application - G ulfstream Quiet Supersonic Jet 



• Conceptual design by Gulfstream Aerospace Corporation (GAC) 

• Mach 1.6 cruise at 45,000 ft 

• Extensive shaping to minimize sonic boom for flight over land 
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Gulfstream Inlet / Rolls-Royce fan 




Gearbox Fairing 


erna 


ompression 


Flow-Through Struts 


Bypass Flow- 

Core Flow s' A 
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Nozzle Shroud 


> 


Bypass Exhaust 
Core Exhaust 


- 0.0 
Mach No. 


Mach 1. 


Rolls-Royce Engine 


Nozzle 


• Isentropic external compression inlet 

• Rolls-Royce engine 

• Plug nozzle 

• Novel bypass duct routes flow internally around large gearbox 
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Computational grids 




Inlet / Fan for a Quiet Supersonic Aircraft 


R. Chima, NASA GRC, Sept. 2009 8 




































































































































































































































































































































































































































































































































































































Axisymmetric inlet / bypass grid 



SIGG grid code by D. Tweedt 

•Algebraic grid with elliptic smoother used for inlet, bypass, and external flow 
• y+ = 1 _ 2.5 


Region 

Type 

Size (x, r) 

Total points 

inlet 

H 

419 x 95 

39,805 

bypass 

H 

439 x 55 

24,145 

external flow 

H 

360x128 

46,080 

Total 



110,030 
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Gearbox fairing blockage 




o.o 

% 0 blockage 


Gearbox fairing modeled using specified blockage in bypass duct grid 
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3-D fan grid 

TCGRID grid code by R. Chima 
•Algebraic H-grid upstream 

• Elliptic C-grids around blades 
•Algebraic O-grids in tip clearance region 

• y+ = 1 _ 2.5 


Region 

Type 

Size (x, q, r) 

Total points 

upstream 

H 

45 x 30 x 95 

128,250 

rotor 

C 

257 x 46 x 95 

1,123,090 

rotor tip 

0 

201 x 15x15 

45,225 

Total 



1,296,565 
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1. Inlet centerbody radius matched to rotor hub. 

2. Tip clearance is unknown. Used 0.025 inches. 
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Grid details 



Blade to blade grid at tip 
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Inlet Corrected Mass Flow [Ib/sec] 


Solution details 





• Solutions run 10,000 - 20,000 iterations at CFL = 2.5 

• Fan pressure ratio and mass flow monitored for convergence 
•7 hr. /10,000 iterations on dual core Xeon CPUs at 3.8 GHz 
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Static pressure c ontours 




• Shocks from fan run upstream into diffuser 
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Mach contours 



•Thick diffuser hub boundary layer —»fan hub radial distortion 
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Core Recovery Bypass Recovery 


Bypass Capture Ratio and Core Recovery 




• Design point: M = 1.7 (over wing) at 45,000 ft. 

• Spillage =1.6 % 

• Bypass recovery = 0.885 

• Core recovery = 0.96 
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Fan Pressure Ratio Adiabatic Efficiency 


Fan Efficiency a nd Pressure Ratio 
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• Hub distortion reduces maximum mass flow by 2.4 % 
and peak efficiency by 0.8 points 
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Percent Span 


Spanwise prof il es at AIP 




•Minimal tip distortion due to bypass design 

• Large hub distortion - uncommon for fans 

• Ps « constant, good BC for isolated inlet calculations 
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Percent Span 



Spanwise prof il es at fan exit 



• Hub distortion persists through the fan. 
Effects on core compressor are unknown. 
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Altitude (k ft) 


Nominal Flight Profile 




Coupled inlet / fan calculations made at each point shown on flight profile 
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8, 0, ft/Qref 


Fan Corrected Conditions 



Mach 
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0 — T 0 1 T q sls 
8 = P 0 / P 0SLS 

m . = rhy/O / 8 


N c = n/Jo 

PR — P ip 

1 lx r 02 ' r O,AIP 


Corrected flow 
Corrected speed 
Pressure ratio 
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Boundary Conditions 

Upstream 

•Subsonic: P 0 , T 0 specified,/? - extrapolated from interior 

• Supersonic: all quantities specified 

Downstream and bypass exit 

• Bypass closed for M < 0.3 

• Subsonic: Ps specified 

• Supersonic: all quantities extrapolated 

Freestream 

• Characteristic perturbation about freestream 

• Reynolds number varies with altitude: Re = f(P 0oo ,T 0oo ) 

Fan 

• Corrected speed varies with altitude: Nc = Q.\fo 

• Q = constant over most of flight profile 

• Q reduced for M > 1.4 due to T3 limit based on cycle deck data from GAC 

• Ps = linear with Nc, matched to desired mass flow at 2 operating points 
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Mach Contours 



The following slides show Mach contours at most points on the flight profile. 
• Note boundary layers at fan face, bypass choking, and shock structure. 
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Mach Contours 
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Mach Contours 







/ 


Takeoff 

• Bypass closed 

• Severe tip distortion 
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Mach Contours 



M = 0.3, 5,000 ft. 

• Bypass open 

• Severe tip distortion 
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Mach Contours 



M = 0.45, 10,000 ft. 

• Diffuser choked 

• Moderate tip distortion 
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Mach Contours 



M = 0.72, A = 18,000 FT 



• Bypass nearly choked 

• Low tip distortion 
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Mach Contours 



M = 0.90, 30,000 ft. 

• Bypass choked 

• Shock on shoulder 

• Low tip distortion 
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Mach Contours 



M = 1.00, 30,000 ft. 


• Possible shock on nacelle 

• Shock on shoulder 

• Low tip distortion 
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Mach Contours 


M = 1.10, 30,000 ft. 

• Shocks on spike, cowl, shoulder 

• Low tip distortion 
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Mach Contours 



M = 1.20, 33,000 ft. 

• Shocks on spike, cowl, shoulder 


• Low tip distortion 
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Mach Contours 



M = 1.40, 39,000 ft. 

• Shocks on spike, cowl, shoulder 

• Low tip distortion 
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Mach Contours 



Cruise, M = 1.60, 45,000 ft. 

• Shocks on spike, cowl, shoulder 

• Hub distortion increasing 


Inlet / Fan for a Quiet Supersonic Aircraft 



R. Chima, NASA GRC, Sept. 2009 35 
































































Mach Contours 



Over wing, M = 1.70, 45,000 ft. 

• Shocks on spike, cowl, shoulder 

• High hub distortion 
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Capture Streamlines 

The following slides show capture streamlines at some points on the flight profile. 

• Note the streamline shape around the bypass splitter lip and the effect on the 
core tip boundary layer. 
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Capture Stream lines 
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Capture Streamlines 




Takeoff 

• Bypass closed 

• Forward motion compresses bypass bubble 

• Streamlines separate at bypass lip 

• Severe tip distortion 
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Capture 


Streamlines 




■ M = OA 5,000 ft. 


• Bypass open 

• High negative incidence on bypass lip 

• Streamlines separate at bypass lip 

• Severe tip distortion 
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Capture 


Streamlines 



M = 1.0, 30,000 ft. 

• Minimum capture area 
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Capture Streamlines 
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Inlet Core Recovery 


Inlet Core Recovery 





Mach 

• Inlet recovery is low at M < 0.7 due to bypass separation 

• Recovery is excellent for 0.7 < M < 1.4 

• Recovery decreases with M > 1.4 due to shock losses 
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Fan Pressure Ratio 


Fan Pressure Ratio 






Mach 

• Fan pressure ratio varies with corrected speed 

• Decreases for M > 1.4 due to reduced Q 
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Fan Pressure Ratio Adiabatic Efficiency 


Fan Pressure Ratio and Efficiency Maps 






•3 speed lines shown for fan with cruise 
distortion profile. 

• Fan operating line is nearly linear except for 
low Mach numbers. 

•Nc varies between 92 and 107 % of design. 

•M = 0.15 - 0.3 is close to the stall line due to 
high tip radial distortion. 

• Fan PR drops with Nc for M > 1.4. 

• Fan operates near peak r| at cruise. 
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Conclusions 



Two coupled CFD codes were used to model the interaction between a 
supersonic Gulfstream inlet and a Rolls-Royce fan. 

Advantages of coupled analysis 

• Eliminates modeling at interface 

• Can track inlet distortion effects all the way through the fan 

• Hard to model isolated components for some cases, e.g. low Mach nos. 

Disadvantages 

• Increased grid / solver complexity 

• Longer solutions times 

Gulfstream inlet / Rolls-Royce fan 

• Design point: good inlet recovery, fan operates near peak efficiency. 

• Hub radial distortion passes through fan, may cause problems for the core. 

• Flight profile: stable operation at all operating points. 

• Identified reduced stall margin at low Mach nos. due to bypass lip separation. 
Lip redesign is being studied. 

•2 papers at AIAA 48 th Aerospace Sciences Meeting, Jan. 4-7 2010, Orlando, FL 
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